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ABSTRACT 

We consider the effect of galaxy intrinsic alignments (IAs) on dark energy constraints 
from weak gravitational lensing. We summarise the latest version of the linear align- 
ment model of IAs, following the brief note of Hirata & Seljak (2010) and further 
interpretation in Laszlo et al. (2011). We show the cosmological bias on the dark en- 
ergy equation of state parameters wo and w a that would occur if IAs were ignored. 
We find that wq and w a are both catastrophically biased, by an absolute value of just 
greater than unity under the Fisher matrix approximation. This contrasts with a bias 
several times larger for the earlier IA implementation. Therefore there is no doubt that 
IAs must be taken into account for future Stage III experiments and beyond. We use 
a flexible grid of IA and galaxy bias parameters as used in previous work, and investi- 
gate what would happen if the universe used the latest IA model, but we assumed the 
earlier version. We find that despite the large difference between the two IA models, 
the grid flexibility is sufficient to remove cosmological bias and recover the correct dark 
energy equation of state. In an appendix, we compare observed shear power spectra 
to those from a popular previous implementation and explain the differences. 

Key words: cosmology: observations - gravitational lensing - large-scale structure 
of Universe - galaxies: evolution 



1 INTRODUCTION 

The gravitational lensing of distant galaxy images has the 
potential to be a powerful cosmological tool. The lensing 
effect directly probes the matter distribution as a function 
of redshift, and thus tells us about the expansion history and 
growth of structure in the Universe. In this way it allows us 
to constrain the dark energy equation of state, or whatever 
is causing the apparent accelerated expansion. 

Weak gravitational lensing poses a number of tough 
technical challenges if its true potential is to be ex- 
ploited. The typical cosmic shear induced on galaxies 
of interest is of order 1%, which is significantly smaller 
than the intrinsic ellipticity of the galaxies themselves. 
Of course we, as observers, have no access to the un- 
lensed galaxy images so we must treat a population of 
galaxies statistically to recover the cosmological informa- 
tion contained in the cosmic shear signal. The correla- 
tion function of galaxy shapes as a measure of gr avita- 
tional lensing w as first prop o sed b y iKaiserl lfl99S ) an d 
first observed by iBacon et al .1 diooch: iKaiser et all ||2000|) : 
IWittman et~all (|2000h and Ivan Waerbeke et all (|200Ch . For 



review s see iBartelmann fc Schneider! (l200ll); iMunshi et ail 



<|2006h : lRefregierl (|2003h . and iHoekstra fe Jairil "(f2008) 



A naive approach to cosmic shear assumes that the in- 
trinsic distribution of galaxy ellipticities is random across 
the sky. If this was the case observed ellipticities on a cer- 
tain patch of sky could be averaged to recover the cos- 
mic shear, because the intrinsic ellipticity would average to 
zero. However it was soon pointed out that this assump- 
tion of random intrinsic ellipticity distribution is unjustified 
dHeayens et al J 120001: ICatelan et alll200ll ; ICrittenden et all 
120011 : ICroft fc Metzlerll2000D . 

In fact galaxies may be expected to align with the large 
scale gravitational potentials in which they form so we ex- 
pect physically close galaxies to be preferentially aligned 
with each other (known as the Intrinsic-Intrinsic (II) cor- 
relation). iHirata fc Seliakl l|2004l ) noted an additional neg- 
ative correlation between foreground galaxies shaped by 
a particular gravitational potential and background galax- 
ies which are lensed by the same potential (known as the 
Gravitational-Intrinsic (GI) correlation) which can be of 
greater magnitude than the II term. 

After the alignment of galaxy ellipticities from linear 
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response to the gravita t ional potential was proposed as an 
effect bv lCatelan et all l|200ll ) the study was put on a firm 
analytic footin g by the introduction o f the Linear Alignment 
(LA) model bv lHirata fc Seliakl (|2004l ), hereafter HS04. This 
approach, in which the orientation of galaxies responds lin- 
early to the large-scale gravitational potential in which they 
form, has become the standard model when including the 
effects of Instrinsic Alignments (IA) . 

Correlated shears of galaxies have been observed at low 
redshift, where they can be attributed to intrinsic align- 
ments, by Brown et al. 2002, and hav e been measured in 
galaxi es s elected to be ph ysical ly c lose in | Mandelbaum et aP 



LVS1C 

( 20061); iHirata et all d2007h: lOkumura et all (|2009h : 



iBrainerd et all (|2009h and lMandelbaum et alj l|2006h . Num- 
ber density - shear correlations are easier to observe since 
the random galaxy ellipticity only ent ers the calculation 
once, an d this has be en cons tra ined in iMandelbaum et al.l 
il2006h; IHirata et all <|2007| ); lOkumura fc Jind i|2009l ); 
IMandelbaum et alj ( 20091 ) and Ijoachimi' et all (|201ll ). 



Hirata et alj (|2007l ) noted that the scale dependence 



of the signal they measured was better matched to theory 
if the non-linear matter power spectrum was used in the 
LA instead of the lin ear power spectrum implied in HS04. 
iBridle fc King! l|2007l ) then used the non-linear matter power 
spectrum in their cosmological forecast calculations. This 
approa ch has been called the Non -Linear Alignment (NLA) 
ansatz. ISchneider fc Bridle! (|2010h attempted a more moti- 
vated solution for assigning power to the IA model at small 
scales through the use of the halo model of galaxy cluster- 
ing. By design their halo model reproduced the LA model 
at large scales. 

The II contribution to intrinsic alignments can be 
taken into account in cosmological analyses by removing 
gala xies with small physical s epara ti on from the analysis 
e.g. iKing fc Schneider! (l2002l. l2003l ); iHevmans fc Heavens! 
(2003); iTakada fc White! (|2004l ). This was attempted in a 
real analysis of COSMOS data in ISchrabback et all (|2009l ) 
by removing the autocorrelation tomograph ic bin . An ex ten- 
sion of this method for GI was suggested in iKind (120051) and 
developed t o a sophisticated level in ljoachimi fc Schneiderl 
(2008, 2009). Alternatively a model may be assumed for all 
the intrinsic alignment contributions, and the free parame- 
ters can be margi nalised over, as demon st rated at th e Fishe r 
matrix level in IBridle fc Kind d2007l); iBernsteinl J2009I); 



Ijoachimi fc Bridle! (|2009l ); ILaszlo et alj (1201 ll); iKirk et alj 



l|201ll ) and demonstrated on real data in lKirk et al, I (|2010l ). 

The redshift evolution of the IA contributions in the LA 
model was found to be incorrect due to a mistake in HS04 
in the conversion between the primordial potential and the 
matter power spect rum. This was correcte d in a new version 
of HS04, issued as IHirata fc Seliakl (|2010l ). hereafter HS10. 
In addition there is an ambiguity in HS04 as to which cos- 
mological epoch is responsible for the "imprinting" of galaxy 
intrinsic alignments. Most starkly the question is, are IAs 
frozen in at some redshift of formation or do they evolve 
with the growth of st ructure, parti c ularly nonlinear clus- 
tering on small scales? iBlazek et al.l (|201ll ) note this as an 
issue in the HS04 approach and estimate i ts effe c t on the GI 
ampli tude as of order 20%. ILaszlo et all (|201ll ). | Kirk et al.l 
l|201ll ) implement one physically motivated solution to this 
questi on (as well as includin g the redshift evolution correc- 
tion of lHirata fc Seliakl feoiO) which we adopt in this work. 



The result of this evolution in the treatment of IAs 
is that much useful theoretical and observational work has 
been conducted using an incorrect implementation of the LA 
model, using often unjustified treatments of the small scale 
seeding of galaxy IAs. In this work our aim is to present 
basic results for the most up to date implementation of the 
LA model for IAs. We present angular power spectra for 
components of the shear-shear (ee), position-position (nn) 
and position-shear (ne) observables as well as the reduction 
in constraining power in measuring dark energy caused by 
a robust treatment of IAs and the biasing of cosmological 
parameters that results from ignoring IAs or treating them 
using an old model. Many of these results reproduce previ- 
ous work in the literature which was conducted using an old 
implement ation of the LA model. S pecifically we reproduce 
Fig. 1 from Ijoachimi fc Bridle! l|2009l ) for the new implemen- 
tation of the LA model which we hope will act as a reference 
for those wishing to apply it in the future. Furthermore we 
investigate the ability of a flexible parameterisation of in- 
trinsic alignments to compensate for using the old LA model 
compared to the latest version we summarise here. 

This paper is organised as follows: In section [2] we de- 
scribe the most up to date implementation of the LA model 
of IAs, detailing our formalism for shear-shear (ee), position- 
position (nn) and position-shear (ne) correlations. Section 
[3] presents the bias on cosmological parameter estimation 
caused by mistreatment of IAs for the latest model, the old 
model and an intermediate model; we then conclude in sec- 
tion 2] In an appendix we summarise the history of the LA 
model and present the major differences between the latest 
implementation and the previous widely used version. 



2 THE LINEAR ALIGNMENT MODEL 

In this section we summarise the latest linear alignment 
model interpretation and its contribution to (i) shear- 
shear correlations (ii) position-position correlation and (iii) 
position-shear cross correlation. 



2.1 Shear-shear correlations, ee 

The Linear Alignment (L A) model for g a laxy i ntrinsic ellip- 
ticities was sugg ested bv | Catelan et al.l l|200ll ) and applied 
quantitatively bv lHirata fc Seliakl 1 20041 ), see section|4]below 
for a more complete history. It assumes that, under gravita- 
tional collapse, dark matter forms into triaxial halos whose 
major axis is aligned with the maximum curvature of the 
large scale gravitational potential. Elliptical galaxies are ex- 
pected to trace the ellipticity of their parent halos, meaning 
that the ellipticity distribution of a population of elliptical 
galaxies will be linearly related to the curvature of the grav- 
itational potential, <j>, via 



c{dl-dl)<p 

2Cdyd X (j>, 



(1) 

(2) 



where C is an unknown constant, fixed for our formalism by 
a normalisation given below. 

In Weak Gravitational Lensing (WGL) we are inter- 
ested in IAs as systematic effects which contaminate a mea- 
sured cosmic shear signal, specifically the shear-shear power 
spectrum. The LA model motivates the amplitude of two 
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separate terms which affect our ability to accurately measure 
cosmic shear. The first, Intrinsic-Intrinsic (II), term arises 
when galaxies are physically close. These galaxies form in 
the same gravitational potential, hence their intrinsic ellip- 
ticities are preferentially aligned. This produces a spurious 
positive correlation which adds to the measured cosmic shear 
signal. 

The second, Gravitational-Intrinsic (GI) , term becomes 
important when treating galaxies which are close on the sky 
but separated in redshift. In this case foreground galaxies 
will align to a foreground gravitational potential which will 
itself contribute to the gravitational lensing of the back- 
ground galaxies. This produces an anti-correlation in ob- 
served ellipticities and subtracts from the expected cosmic 
shear signal. 

In terms of projected angular power spectra we can 
write the measured ellipticity power spectrum, C" , as a sum 
of the gravitational lensing shear power spectrum and two 
IA terms 



C"{l) = C Ga {l)+C 11 (l) + C ul (0 



GI i 



(3) 



The cosmic shear power spectrum, Cy (I), under the 
Limber approximation, is written as the integrated product 
of the matter power sp ectrum, Pssjk, y), and the lensing 
weight function, Wi(x) (|Hulll999l '). 



C GG (l) 



X 2 



(4) 



where Pss(k, x) = Pss ( f K \ x ) i x)> X i s the comoving distance 
along the line of sight, }k{x) is the comoving angular diame- 
ter distance, i,j denote a tomographic redshift bin pair and 



3 H () Q, m X 



2 c 2 a 



dx'Mx 



i^X -x 



(5) 



with Ho the Hubble parameter today, Q m the dimensionless 
matter energy density, c the speed of light, a the dimension- 
less scale factor and rii(x') is the galaxy redshift distribution 
for a particular bin i. 

Similarly, the projected angula r power spectra for the 
I A terms are (|Hirata fe Seliakll2004 ) 

O) = ^Mx)nj( X )Pu(k,x) (6) 

Jo X 

C?/(l) = f Xh °' ' %W i {x)n j { X )PGi{k,x), (7) 
Jo X 

where, unsurprisingly, the II correlation depends only on the 
galaxy redshift distribution, n(x), while GI depends on the 
product of the galaxy redshift distribution and the lensing 
weight function. 

The IA "power spectra", Pn(k, x) an d PGi{k,x), are 
unknown functions of scale and cosmic epoch. In the case of 
the LA model they become 

Pu(k, X ) = (-C7 lP (x = 0)) 2 P^ n (fc, X = 0) (8) 

P G i(k, X ) = -C lP ( X = 0)yJP$(k, X = 0)^|i^(9) 

where Pss is the matter power spectrum, p(x = 0) is the 
matter density today and D(z) is the linear growth factor as 
a function of redshift z, where it is implicitly meant that z = 
z(x)- Ci is the normalisation of the I A contribution and we 



14 (/i 2 M Mpc- 3 )- 1 , following 
who match to the power spectra in 



use a fiducial value o f 5 x 10 
Bridle fc Kingl J2007h who r____ 

Hirata fc Seliakl (|2004) which are based on SuperCOSMOS 



Brown et all (|2002l ). 



data from 

The II term is related to the linear matter power spec- 
trum, PsTj because we assume that galaxy intrinsic ellip- 
ticity is imprinted at the (early) epoch of galaxy formation 
and subsequently "frozen in". As such the power spectrum 
of its distribution does not undergo the nonlinear evolution 
that the distribution of matter clustering does at late times. 
The factor of \/Pgs L {k,x = 0)Pss(k,x)/D(z) in the GI term 
reflects the dependence on both the intrinsic ellipticity dis- 
tribution and the full, late time, mass distribution via grav- 
itational lensing. The exact relation between intrinsic ellip- 
ticity, the epoch of galaxy formation and the evolution of 
galaxy clustering is a vexed and disputed one. A more de- 
tailed discussion of the history of the treatment of this re- 
lationship can be found in the appendix below. We believe 
that our approach best reflects the physical understa nding 
of the LA model as proposed bv lHirata fc Se liak ( 2004) who 
related the intrinsic ellipticity distribution to the Newtonian 
potential at the time of galaxy formation. 



2.2 Position-position correlations, nn 

Any cosmic shear survey provides not only a catalogue of 
measured galaxy ellipticities but also a record of galaxy po- 
sitions through their location on the sky and an estimate of 
their redshift. Analogously to the cosmic shear power spec- 
trum we may define the galaxy position density power spec- 
trum, C""(l), the Fourier transform of the galaxy position 
two-point correlation function. It is conceptually useful to 
divide the contributions up as follows, 



= c"(0 + c mm (i) + c 9m (i). 



(10) 



Here we have separated the contribution from galaxy clus- 
tering itself, gg, from the contribution due to lensing mag- 
nification, mm, and the cross-correlation of the two, gm. 
The lensing contribution manifests because a galaxy may 
increase (or decrease) in size as a result of the image distor- 
tion. The application of Liouville's theorem tells us that the 
process of WGL conserves surface brightness density hence 
a larger (smaller) galaxy will appear brighter (dimmer) due 
to the lensing effect. For a magnitude limited survey this 
can affect the statistics of galaxy clustering as otherwise too 
faint galaxies are promoted into the survey by magnifica- 
tion and vice-versa. As the nn correlation is effectively just 
a counting of galaxies and their relative positions there is no 
ellipticity contribution and we can ignore IAs. 

The equation for the pure galaxy clustering term is rel- 
atively straightforward in the Limber approximation 

- d X 



X 



L ni(x)nj(x)b g (k, z)P ss (k,x), (H) 



where the window function is the galaxy redshift distribu- 
tion as we would expect and we have introduced the galaxy 
bias term, b g , to reflect our understanding that galaxies are 
a biased tracer of the underlying matter distribution. As 
before, Pss(k,x) = Pss( f K \ x ) > x)- In general b g is an un- 
known function of scale and cosmic epoch. The details of 
the galaxy bias formalism we assume in this paper, along 
with other bias terms, are explained in section [3] below. 
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The magnification power spectrum and galaxy 
clustering-magnification cross-term can be written 

arm = ^-lx^-i^f « (12) 

Cf™(l) = 2(a J -l)C? G (I), (13) 

where aj is defined as the slope of the luminosity function, 
evaluated at the median redshift of the bin i, and C gG (l) 
is defined in sectio n J2.31 In this approach we are following 
Ijoachimi fc Bri dle (2009). We treat each of the a.i terms as 
a free parameter in the model, an d take fiducial parameters 
from eqns. 32, 33 and Table 2 of lJoachimi fc Bridle! (|2009l ) 
(see Appendix A of that paper for more details). 

2.3 Position-shear correlations, ne 

As well as ee and nn correlations themselves, we can cross- 
correlate the two fields to form the galaxy positio n-shear, ne, 
cross- correlation functions as first proposed by iHu fc Jam! 
|2004h in the context of dark energy. This is often referred to 
as galaxy-galaxy lensing, in which the mass of a foreground 
galaxy distorts the shape of a background galaxy. Here we 
consider the general cross-correlation which includes contri- 
butions from larger dark matter structures, from magnifi- 
cation and also from intrinsic alignments if the correlated 
galaxies are physically close. We can write the contributions 
to the full ne term as 

C nt (l) = C gG {l) + C 9l {l) + C mG (l) + C mI (l) (14) 

|joachimi fc Bridle! [2009^ . 

The individual expressions for the terms in the C ne (l) 
expansion contain combinations of quantities already con- 
sidered 

cfd) = f hor ^MxWx) 

Jo X 

bg(k,z)r g (k,z)Pss(k,x) (15) 
Jo X 2 

b g (k, z)r g (k, z)bi(k, z)n(k, z)P G i{k, x) (16) 

C% G (l) = 2( Qi - l)Cg G (Z) (17) 

Of(/) = 2(^-1)^(0, (18) 

where we have introduced the additional free functions r g 
and n which appear due to the possible stochastic re- 
lation between the g alaxy and dark matter distributions 
|Dekel fc Lahavl 1 19991 ) as discussed in section [3] below. As 
before Pss{k,x) = Pss{ f K \ x ) > X)- We note that Poi(k,x) 
appears in the C?£ (Z) term because it describes how the 
dark matter distribution relates to the intrinsic alignments, 
and we assume the relation between the galaxy distribution 
and the dark matter distribution can be accounted for by 
the product b g r g . 

2.4 Summary of Observables and Fields 

The observables considered in this paper and the different 
fields which contribute to each observable are summarised 
in Table Q] In Fig. \T\ we plot the angular power spectra of 
all considered components for a fiducial stage-IV lensing sur- 
vey. For our fiducial model in this paper we use the following 



Observables 




cc 


shear-shear 


till 


position-position 


TIC 


position-shear 


Fields 




G 


gravitational lensing 


I 


intrinsic alignment 


9 


galaxy clustering 


m 


cosmic magnification 



Table 1. Summary of observables considered in this paper and 
different fields which contribute to each observable. 

survey specification: 20,000 deg 2 with a galaxy number den- 
sity of 35 arcmin -2 and a redshift distribution given by the 
Smail-type n(z), 

n(z) = z a eiq ? (- (£) J, (19) 

with a = 2, ft = 1.5, zq — 0.9/\/2, divided into 10 tomo- 
graphic bins with equal number density out to redshift 3 
and normalised so that J n(z)dz = 1. A Gaussian photo- 
metric redshift error of a z — 0.05(1 + z) is assumed with no 
catastrophic outliers in redshift. We constrain the set of cos- 
mological parameters {Q, m , wo,w a , h, as, Qb, n a } which take 
the values £l m = 0.25, Wo = — 1, W a = 0, h = 0.7, as = 0.8, 
Qb = 0.05, n s — 1. All results assume a flat, ACDM cosmol- 
ogy. Linear matter pow er spectrum is calculated using the 
lEisenstein fc Hrj (|l998T ) fitting function and n onlinear cor- 
rectio ns are applied where appropriate using the lSmith et al.l 
(|2003h formalism. 

Fig. [J is presented and form atted in the same way as 
Fig. 3 in lJoachimi fc Bridle! ([2009). The purpose is to supply 
a new reference plot using the most up to date implementa- 
tion of the LA model for I As. All lines which do not include 
an "I" ter m will be identical in both versions of the plot. 
However, as Ijoachimi fc Firi dlc (2009) use the original HS04 
implementation, f rom before the publica tion of the erratum 
incorporated into lHirata fc Seliakl (|2010T ). and present lines 
for the NLA ansatz, there are differences in any lines con- 
taining an I contribution. The NLA version of HS04 is the 
most widely used in the literature to date; we refer to it 
subsequently as HS04NL. Note that that we are employing 
the NLA prescription in HS04NL e ven though this was not 
explicitly discussed in the original iHirata fc Seliakl (|2004T ) 
paper. See the appendix for a more detailed explanation of 
the evolution of the LA model and the differences between 
various implementations. 

The upper triangle of Fig. [T] shows the contributions to 
the shear-shear correlation function in black for every other 
tomographic bin pairing. As usual the lensing contribution 
is the largest on all angular scales, but is most dominant 
in the auto-correlations at high redshift. The II intrinsic 
alignment term is most important in the auto-correlations 
and negligible for widely spaced bins. The GI contribution is 
largest for separated bins (shown as IG for compactness - see 
caption) . The intrinsic alignment contributions are typically 
one to ten per cent of the lensing signal. 

The contributions to the position-position correlation 
function are shown in pink in the upper triangle and as 
expected, the auto-correlations are dominated by the in- 
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100 1000 10 100 1000 10 100 1000 10 100 1000 

till 



100 1000 

I 



Figure 1. Fiducial power spectra for all correlations considered. The upper right panels depict the contributions to the ee (in black) 
and nn (in magenta) correlations. The lower left panels show the contributions to correlations between number density fluctuations and 
ellipticity. Since we only show correlations C^^(l) with i ^ j, we make in this plot a distinction between ne (in red; number density 
contribution in the foreground, e.g. gG) and en (in blue; number density contribution in the background, e.g. Gg) correlations. In each 
sub-panel a different tomographic redshift bin correlation is shown. For concision only odd bins are displayed. In the upper right panels 
the usual cosmic shear signal (GG) is shown as black solid lines; the intrinsic alignment GI term is shown by the black dashed lines; 
the intrinsic alignment II term is shown by the dotted black lines; the usual galaxy clustering signal (gg) is shown by the magenta solid 
lines; the cross correlation between galaxy clustering and lensing magnification (gm) is shown by the magenta dashed lines; the lensing 
magnification correlation functions (mm) are shown by the magenta dotted lines. In the lower left panels the solid blue lines show the 
correlation between lensing shear and galaxy clustering (Gg); the blue dashed lines show the correlation between lensing shear and lensing 
magnification (gm); the blue dot-dashed lines show the correlation between intrinsic alignment and galaxy clustering (Ig or equivalently 
gl); the red solid lines show the correlation between galaxy clustering and lensing shear (gG), which is equivalent to the blue solid lines 
with redshift bin indices i and j reversed; similarly the red dashed lines show the correlation between lensing magnification and lensing 
shear (mG), for cases where the magnification occurs at lower redshift than the shear (i < j); finally the dotted lines show the correlation 
between lensing magnification and intrinsic alignment (ml). 
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trinsic clustering of galaxies. The cross-correlations of sep- 
arated bins can be dominated by the cross-correlation 
of galaxy clustering and magnification. As expected, the 
magnification-only term is most important at high redshift, 
although it never exceeds the other terms for the redshift 
ranges we consider. 

The lower triangle shows the various contributions to 
the position-shear cross-correlation function, and shows the 
effect of reversing the tomographic bin order where relevant 
(see caption). The correlation between lensing shear and 
galaxy clustering is the strongest of these contributions and 
is largest when the tomographic bins are separated with the 
galaxies in the foreground. This is the usual galaxy-galaxy 
lensing contribution. Interestingly it is still large in the auto- 
correlations which is due to the galaxies in the nearest part 
of the redshift bin lensing the galaxies in the farthest part. 
There is a non-negligible contribution in some of the cross 
terms from the lensing of galaxies in a nearer tomographic 
bin by galaxies in a farther tomographic bin, due to the 
overlap in redshift distributions (solid blue lines). At high- 
est redshifts the magnification terms start to dominate. The 
cross-correlation between magnification and intrinsic align- 
ment is mostly one of the smallest terms. 



3 BIASING OF COSMOLOGICAL 
PARAMETERS 

The LA model is currently our best description of galaxy 
intrinsic alignment over a range of cosmologies. It has been 
shown that ignoring IAs in an analysis can significantly bias 
the measurement of fundamental cosmological parameters 
jBridle fc Kind 120071 ). In this section we quantify this bias 
and examine the effect of moving from the old standard LA 
implementation to the new one. We also explore the bias 
produced by using the old implementation rather than the 
new. 

A naive approach to cosmic shear measurements would 
ignore IAs and measure values for cosmological parameters 
which are systematically biased. By employing a flexible IA 
model and marginalising over a set of nuisance parameters 
we lose precision but hope to produce unbiased cosmological 
measurements. In this section we quantify the cosmological 
bias which results from ignoring IAs and also the cosmo- 
logical bias which results from employing the HS04NL LA 
model rather than the current LA model. 

To this end we em ploy the c osmological bias formalism 
of iHuterer et al. | (I2006T) (see alsolAmara fc Refregieil {2006) 
and appendix A of lJoachimi et al.l 1)20111 )). 



= ]T Adj (Cov [£7 y (0, O^G)])" 1 



(20) 

where 5p a is the cosmological bias on each parameter consid- 
ered, dj are the projected angular power spectra, dC ^^ 
are the derivatives of these power spectra with respect to 
the cosmological parameters, Cov [dj(l), C mn (0] is the co- 
variance mat rix of the power spectra, calculated according 
to eqn. 39 in IJoachimi fc Bridle! (2009), and F a p is the in- 
verse of the Fisher Matrix (FM) for our set of cosmological 



parameters, calculated as 

Each of the FM terms is calculated using the assumed IA 
model. Adj is the difference between the data vector cal- 
culated using the assumed model and the true model, 



ACV 



^assumed 



, '(*)• (21) 
Fig. [2] shows the cosmological bias on the dark energy 
equation of state parameters if IAs are assumed not to exist 
but they are truly present according to one of: the HS04NL 
implementation; the new model; or an intermediate between 
the two, called HS10NL. The bias is the distance between 
the centre of the "true" contour [-1,0] and the centre of the 
contours which assume no IAs. 

Note that the bias parameterisation used is based on 
the FM formalism which assumes purely Gaussian errors in 
the parameters and observables which respond linearly with 
respect to the parameters in question. The parameterisation 
is very likely to break down for parameter values more than 
~ 2 — 3cr away from the fiducial cosmology. As such the 
more extreme bias values seen in Figs. [2] and [3] should not 
be read as exact. What they can tell us is rough relative bias 
between different scenarios and it is clear that any scenario 
producing an absolute bias on wo of order unity or above can 
be considered to be "catastrophically biased". For clarity 
we use the notation ctfm when quoting the size of biases 
on parameters with respect to the errors on the unbiased 
constraints to remind the reader that they are calculated 
assuming the Fisher Matrix formalism is valid. 

The intermediate model (HS10NL) is introduced to dis- 
entangle the two effects which change between HS04NL and 
the new LA implement ation. It includes the correct fac- 
tors of a, introduced by iHirata fc Seli aid l|2010h and shown 
in equations 1381 but it always applies the non-linear mat- 
ter power spectrum, Ph n (k,z), to the IA terms, following 
the lBridle fc Kind (|2007f l interpretation of HS04NL. So this 
HS10NL model applies the correct redshift evolution to IAs 
but assumes that the effects of nonlinear clustering are al- 
ways present, even in the II term, by adopting the NLA 
ansatz. 

Ignoring IAs causes strong biasing on the dark energy 
parameters, no matter which IA model is assumed to be 
true. The new implementation is the least biased at ~ 8<tfm 
away from the true model. However, ignoring IAs appears 
to still bias wo by of order ~ 1.5. The HS10NL model is 
more biased at ~ 20<tfm- The HS04NL implementation pro- 
duces the strongest cosmological bias, with the contour far 
outside the plotted area. At a point this far away from the 
fiducial parameter values the Gaussian assumptions of the 
Fisher matrix and bias formalisms certainly break down and 
it would be wrong to put much faith in the exact direc- 
tion/distance of the predicted bias. What is clear however is 
that the effect is very strong. This finding is what we would 
expect given that the move from HS04NL to the new imple- 
mentation has reduced the impact of IAs, not only through 
changes to their redshift evolution, but also through the re- 
moval of IA power on small scales. 

The general biasing trend is the same when we calcu- 
late for a stage-Ill type survey, like the Dark Energy Sur- 
vey (DES), with biases of ~ 4<tfm for the latest model, 
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Figure 2. 95% confidence limits on the dark energy parame- 
ters wo and w a for the IA implementation from this work [blue 
contour] and the biased constraints if it is assumed that IAs do 
not exist but in fact the new implementation [black contour], the 
HS10NL implementation [green contour] or the HS04NL imple- 
mentation [off plot, direction indicated by red arrow] are true. The 
cosmological bias is given by the offset between each displaced 
contour and the fiducial values of (—1,0) [black cross]. Results 
are presented for an observable data vector of shear-shear (ee) 
correlations only. It is assumed that we enjoy perfect knowledge 
of the IA contribution for each implementation. 



y 

w 

ee, marginalised over IA grid 




EE + n£+nn, marginalised over IA grid 




_ Latest IA model in data 
latest IA model in fit 
Latest IA model in data 

"HS04NL in fit 



Figure 3. 95% confidence limits on the dark energy parameters 
wo and w a for the IA implementation from this work [blue con- 
tours] and the biased constraints if the old HS04NL implementa- 
tion is assumed [red contours] . The cosmological bias is given by 
the offset between the displaced contour and the fiducial values of 
(—1,0) [black dotted lines]. Results are presented for an observ- 
able data vector of shear-shear (ee) correlations alone [left-hand 
panels] and the full combination of shear-shear, shear-position and 
position-position (ee + ne + nn) correlations [right-hand panels]. 
Each probe combination is shown for the case of perfect knowl- 
edge of IAs and galaxy bias [top panels] and for the case where our 
ignorance is accounted for by marginalisation over the nuisance 
parameter grid with N k = N z = 7 [bottom panels]. 



~ 8o"fm for the HS10NL implementation and ~ 30<tfm for 
the HS04NL model. 

So far we have not taken into account any uncertainty 
in the IA model, but assumed they are zero in the param- 
eter fitting, whichever model we employ to describe them 
in the true model. In reality, we are aware that our knowl- 
edge of IAs from simulations and observations is still de- 
veloping and relatively uncertain. The case is similar for 
galaxy bias, introduced in section 12.21 and employed be- 
low. To parameterise our ignorance of both effects and their 
cross-correlations we use a grid of nuisance parameters 



X = A x Qx(k,z), 



(22) 



where Ax is a constant amplitude parameter, free to vary 
about a fiducial value of 1, and Qx(k, z) is a grid of Nk x N z 
nodes logarithmically spaced in k/z space, each of which 
is allowed to vary independently around a fiducial value 
of 1. A final smooth grid is created by spline interpola- 
tion over the values of the grid nodes. For more details o f 
this nuisance parameter grid see iJoachimi fc Bridle! (2009), 
iLaszlo et ail j201ll ), [Kirk et all l|201ll) . This was inspired by 
the mar ginalisation of ? which l ed to the grid implementa- 
tion in IJoachimi fc Bridle! (|2009l V 

We define four sets of nuisance parameters which each 



take the form of this grid, 

bi(k,z) = A bl Q bl (k.z) 

b g {k,z) = A bg Q bg (k.z) 

r r (k,z) = A ri Q ri (k.z) 

r g {k,z) = A r Q r (k.z). 



(23) 
(24) 
(25) 
(26) 



The "6" terms can be understood as bias terms between 
the power spectrum of a field X and the matter power spec- 
trum, i.e.: 



Pxx(k,z) _ , , 



Pss(k, 



b x (k,z) 



(27) 



and "r" terms are biasing terms which arise in cross- 
correlations, i.e.: 



Px Y {k,z) 
Pss{k,z) 



= bx (k, z)r x {k, z)b Y (k, z)r Y {k, z). (28) 



The "6" and "r" terms are both f unctions of the c l usteri ng 
and stochastic bias as described in lDekeffc Lahavl(|l999l ). 

These sets of nuisance parameters then appear, in differ- 
ent combinations, in each angular power spectrum integral 
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where galaxy clustering, g, or the IA term, I, appear: 



(29) 



%n z { X )n ] (x)b]{k,z)P II {k, X ) 
X 2 

%W 1 (x)n 3 ( X )b I (k, z)r r (k, z)P G i(k, X ) (30) 
X 2 

X 2 



}rn{x)nj(x)b g (k, z)P ss {k,x) 



(31) 



f^rn ( X ) W 3 ( X )b 3 {k,z)r g (k,z)P SG {k, x ) (32) 
J X 

/ ^ n i(xX(x)MM)ri(&>«) 



bg{k,z)r g (k, z)Psi{k,x), 



(33) 



while nuisance parameters appear in C^ m (Z) and Clj 1 (I) due 
to their dependence on Cf^ (I) and Cfj 1 (I) respectively. In 
general this means that, for the full suite of observables 
ee + ne + nn, we marginalise over a total of 4(1 + Nk x N z ) 
nuisance parameters, plus 30 that parameterise the magnifi- 
cation effects. As before Pss(k, x) = Pss{ f K \ x ) iX)- When we 
employ the nuisance grid in this work we set Nk = N z — 7, 
giving a total of 230 nuisance parameters . This follows the 
practice used in Ijoachimi fc Bridle! ()2009l ) for the purposes 
of comparison (see their Table 2). We also restrict the in- 
formation used in the galaxy field to linear scales as in 
iRassat et all (|200Sh . 

This current paper presents the most advanced inter- 
pretation of the LA model. However, much work of the 
work on IAs in the literature has been conducted using a 
non-linear version of th e original interpretation presented in 
iHirata fc SeTiakl (|2004l l. In Fig. [3] we quantify the bias intro- 
duced on the dark energy parameters from the use of this 
HS04NL approach rather than the latest implementation in- 
troduced in this work. When the IA signals (whether from 
the old or new models) are assumed to be known perfectly 
there is a very strong biasing effect which would result in a 
very worrying systematic mis-measurement of cosmology. 

In the case where ee correlations alone are considered, 
the bias is ~ 25o"fm of the constraining power of the cosmic 
shear signal given the true I A contribution. When the full 
set of correlations including position information, ee + ne + 
nn, are considered this biasing increases to ~ 65<tfm- As 
previously discussed, when the effect is this strong the exact 
position/direction of the biasing is likely to be well outside 
the competence of the FM-based formalism but the fact that 
there is a very significant effect should not be doubted. 

When we employ the robust parameterisation of both 
IA and galaxy bias uncertainties via the grid of nuisance 
parameters described above we see a striking change in the 
results. As we would expect the constraining power of either 
probe combination decreases as we have marginalised over 
130 free parameters which represent our ignorance of the 
exact details of IAs (and another 100 for galaxy bias) . What 
is also clear is that the biasing of our cosmological estimates 
has decreased to well within the la error of the respective 
probes. 

The constraints on wo and w a from the full ee + ne + nn 
combination after marginalisation over 230 nuisance param- 
eters are only reduced by a factor of two compared to the 
naive (and heavily biased if the wrong IA model is used) 
constraints from ee alone when it is assumed IAs are per- 



fectly known. This corresponds to a factor of 2.6 reduction 
in DETF FoM l|Albrecht et al.ll2006l ). 



4 CONCLUSIONS 

Galaxy IAs are the most important physical systematic in 
the study of cosmic shear. As increasing volumes of data 
become available from WGL surveys more interest is be- 
ing paid to their correct treatment. Much of this work has 
been dominated by the LA model, originally introduced in 
HS04 (and the NLA extension which we refer to as HS04NL). 
This original im p lemen tation was subsequently corrected in 
IHirata fcS cliak (2010) and detailed attention paid to th e 
treatment of non -linear clustering in iLaszlo et al.l l|201ll ); 
iKirk et all (|201ll ). 

However, much of the existing literature has been pro- 
duced using the uncorrected HS04NL model. In this paper 
we have provided a brief explanation of the evolution of, and 
context surrounding, the LA model for IAs, highlighting the 
most important differences between HS04NL and the latest 
implementation. 

We have calculated the angular power spectra of the 
cosmic shear observables , correcting the implicit mistake of 
Ijoachimi fc Bridle! (2009), which was based on HS04, which 
we hope will act as a new reference for those interested in 
applying IAs to their cosmic shear analysis. 

The main motivation for the study of IAs is the mea- 
surement of unbiased constraints of key cosmological pa- 
rameters. We show that the new LA implementation signif- 
icantly reduces the impact of IAs, and hence the bias, but 
that the effect is still very significant, producing a bias at 
the tens of a level when we assume perfect knowledge of IAs. 

If we did know the IA signal perfectly then we could 
produce unbiased measurements by simply subtracting the 
IA signal from our measured cosmic shear. In practice it is 
useful to parameterise our ignorance of the true IA signal 
through a set of nuisance parameters which are marginalised 
over to produce weaker but hopefully unbiased cosmological 
estimates. We show that a robust grid of 130 nuisance pa- 
rameters for IAs and magnification uncertainties, allowed to 
vary in scale and redshift, effectively removes the bias due 
to assuming an incorrect IA model. 

The same effect holds when our observables are ex- 
tended to include shear-shear, position-position and shear- 
position correlations. The extra observables increase con- 
straining power so that we are able to produce unbiased 
constraints on ee + ne + nn which recover 40% of the con- 
straining power of the highly biased constraint from ee alone 
when the wrong IA model is assumed. 

The LA model has allowed a firm foothold on the study 
of IAs to develop over the last decade. We have detailed an 
updated implementation of the most used IA model and its 
physical motivation, showing a reduction in biasing of cos- 
mological constraints. In addition we reiterate that a robust 
nuisance parameter model can control the biasing due to IAs 
for either the old or new implementations. Together, these 
results should give us confidence that the IA effect is un- 
der control as we begin to analyse the first data from large 
cosmic shear surveys. 

However our somewhat pessimistic approach carries the 
cost of reduced constraining power. Future work, focused 
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on accurate simulation and measurement of IAs is sure to 
provide more detail on the physical mechanisms responsi- 
ble for the initial intrinsic ellipticity distribution and its 
evolution. This better knowledge of IAs will improve our 
ability to model them and reduce our dependence on brute- 
force marginalisation over nuisance parameters. As such the 
marginalised constraints which we present in Fig. [3] may be 
a worst-case scenario. With better knowledge of IAs, better 
constraints on cosmology will be possible. 
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APPENDIX: HISTORY OF THE LA MODEL 

The hi story of galaxy intrinsic alignment studi es goes 
back t o |j. Buergers. H. van de Hulst H.I |l95ll ) and ISciamal 
l| 19551 ) who discussed that galaxies may acquire angular 
momentum through tidal torquing, in which anisotropic 
shear flows distort the protogalaxy and lead to a lo- 
cal rotation. The theory was further developed a nalyti- 



cally by iPeebles! (|l969l ). iDoroshkevichl Jl970l) and 
l|l984l ). and first simulated by iHeavens fc Peacock] 



White! 
19881) . 



If galaxy orientations are correlated with their angular 
momenta then we may expect neighbou ri ng galaxies to 
be ali g ned (e.g. van den Bosch et al.ll2002l). IHeavens et al.l 
i|2000l ); ICroft fc Metzlerj |2000l ); ICatelan et all (|200 ll ) simul- 
taneously pointed out that this could lead to a contami- 
nating term in weak gravitational len s ing, a n d this was in- 
vestig at ed fu rt her in Crittenden et al. (2001); Mackcv et all 
2002); Ijind (I2002TI; iHevmans et all il2004h; IHirata et all 



20041 ) and lBridle fc Kind (|2007l ). See ISchaefeij |200Sl) for 



a m ore recent review. 

ICiotti fc Duttall| 19941 ) described a second ansatz to mo- 
tivate the intrinsic alignments of elliptical galaxies, based on 
tidal stretching of the host halo by the surrounding large- 
scale structure. ICatelan et ail (|200ll ) showed that the corre- 
lation function of the intrinsic alignments should be propor- 
tional to that of the tidal field and thus be significant even 
on large scales. This contrasts with the correlation function 
expected from tidal torque theory which has a higher order 
dep endence on the tidal f i eld and thus decays mor e rapidly. 

IHeavens et ail (|200fj| ); ICroft fc Metzlerl |2000l ) used the 
shapes of dark matter halos in A^-body simulations as a 
proxy for elliptical galaxy shapes and estimated the pos- 
sible contamination to the (then) recent cosmic shear de- 
tections. They found the contamination to be a fraction of 
the measured signal but IHeavens et all |2000r i commented 
that low redshift surveys would be strongly affected and 



ICroft fc Metzleil (|2000D point out that future tomographic 
mea surements will a l so be significantly affected. 

IHirata fc Seliakl (|2004T ) (which we are calling HS04) dis- 
cussed both models and used the term 'linear alignment 
mode! to describe the alignment of the galaxy halo with 
the tidal stretch of the local gravitational potential. Previ- 
ous to HS04 the cosmic shear contamination was believed 
to come from physically close galaxies due to both galaxies 
forming in the same large-scale structure. HS04 identified 
a second term due to the intrinsic alignment of one galaxy 
near to a large mass which would correlate with the grav- 
itational alignment of a distant galaxy lensed by the same 
mass. The first of these two was labelled 'intrinsic-intrinsic' 
(II) correlation and the second 'gravitational-intrinsic' (GI) 
correlation. 

Following the approach of ICatelan et all i|200ll ). HS04 
assume that galaxy intrinsic ellipticity follows the linear re- 
lation 

7/ = ~4^ (A '~ A ^' 2A:cAa)5[<I>p] ' (34) 
where <E>p is the Newtonian potential at the time of galaxy 
formation, smoothed by some filter 5 that cuts off fluctu- 
ations on galactic scales, G is the Newtonian gravitational 
constant, A is a comoving derivative and Ci is a normalisa- 
tion constant. 

After relating the primordial potential to the linear 
matter power spectrum, HS04 calculate the II and GI power 
spectra 



D HS04/, n 

Pi i {k,X! 



Pgi (k,x) = 



-Cip(z) 



D(z) 
dp(z) 



Ps5 n (k, X ) 



D(z) 



Pss n (k, X ) 



(35) 



(36) 



in terms of the linear theory matter power spectrum at co- 
moving distance x corresponding to redshift z, P^ 1 5 n (fc,x)- 
Here p(z) is the mean density of the Universe, D(z) = 
(1 + z)D(z) is the rescaled growth factor D(z). 

Note that HS04 actually present the power spectra for 
71, the density weighted intrinsic shear (and 5, ji for the GI 
term). This is because galaxies are not randomly positioned, 
but are expected to form preferentially in regions of higher 
matter density. This leads to higher order terms in the IA 
power spectra which are usually ignored because they are 
about an order of magn itude smaller than th e leading term 
in the above equations (|Bridle fc Kindl2007h . 

These equations were used in several publications un- 
til the discovery of missing factors of a 2 which arose in 
the conversion factor between the density perturbation and 
the primordial potential. T hese were corrected in an erra- 
tum lffirata fc Seliakl (|2010l ). which we call HS10. In the up- 
dated version, eons. [35l and 1361 become 



D HSl(),-7 N 

Pi i {k,X! 



Pgi (k,X) = - 



-Cip(z) 



D(z) 
Cip(z) 



a 2 ) P l sT(k,x) 



D(z) 



a 2 PsT(k,x)- 



(37) 



(38) 



Originally HS04 related their Pu and Pgi to the 
linear matter power spectr u m Pi ". It was subsequently 
indicated by IHirata et all |2007) that more power at 
small scales may fit the data better and IBridle fc Kin3 
(2007) proposed that the non-linear matter power spec- 
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trum be substituted for the linear in what became known 
and the Non-Li near Alignment (NLA ) model, which was 
then a pplied in iBridle fe Kind ([200711: ISchneider fc Bridle! 
d2010h: iMandelbaum et all (120091 ); jjoachimi et all d201ll); 
iKirk et al.1 (|2010l ): iLaszlo et alj |201ll ) and iKirk et all 
1 20111 ). It is this common NLA approach, which we describe 
as HS04NL, 



p£ S04NL (fc,x) 



D HS04NL/, \ 

Pgi (k,X) 



-dp(z 



D(z) 
Cipjz] 
D(z) 



Pss(k, X ) 
Pss(k,x), 



(39) 



(40) 



where the non-linear matter power spectrum is used in both 
equations. 

In this work we employ a more physically motivated 
approach based on the idea that galaxy intrinsic alignments 
are seeded at some early epoch of galaxy formation and do 
not evolve w ith redshift. This wa s the original intent of the 
LA model in lCatelan et al.l (|200ll ) and HS04. This principle 
is apparent from the original HS04 equations by the presence 
of the 1/D(z) terms which divide out the growth. Therefore 
we relate our II term directly to Pl 6 n but the GI term picks 
up some contribution from the later (linear and nonlinear) 
growth of large scale structure via its dependence on cos- 
mic shear. This is included in the formalism by the geomet- 
ric mean of the linear and nonlinear matter power spectra 
in equation 9l Thi s was proposed an d discussed further in 
ILaszlo et alTi|201ll ): IKirk et all l|201ll ). 



Note that equations 1391 and 1401 are the same as our ex- 
pressions in eqns.|8] and [9] with the prefactor in brackets sim- 
plified, and there is a different treatment of the linear/non- 
linear matter power spectra. 

An improvement on the contributions at non-linear 
scales was developed bv lSchneider fc Bridle! (|201p| ) by devel- 
oping a halo model of intrinsic alignments. We do not use it 
here as it does not contain full flexibility of the cosmological 
model. 

In this paper we discuss and illustrate the difference be- 
tween the models, summarised in Fig. [4] The changes can 
be divided int o two types: (i) the cor rect use of the a 2 terms 
as included in lHirata fc Seliakl (|2010l ) and (ii) the consistent 
use of the linear matter power spectrum for the II correla- 
tion but the square root product of the linear and nonlinear 
power spectra for the GI correlation (see Equation [5} . The 
upper panel of Fig. [4] includes only the first change by plot- 
ting the square root of the Ci integrand without the contri- 
bution from the matter power spectrum. The terms plotted 
for bins i = 1, 10 are 



II 



GG 

J-J-HS04 
This Work 



2wGp(z)a 2 
7i 



-Wi( X ) 
Cip{z)am(x) 
Cip(z = Q)m(x)- 



(41) 
(42) 
(43) 



The lower panel includes the square root of the matter power 
spectrum contribution, plotting 



II 



GG 

JJ-HS04NL 
This Work 



C 1 p(z)an i ( X ) ^/Pss(k,x) 



Cip{z = 0)m(xWPS?(k, x ). 



(44) 
(45) 
(46) 



Correct inclusion of the a 2 terms produces an IA signal 
which is constant with redshift, rather than one which in- 
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Figure 4. Integrand of the GG angular power spectrum at a mul- 
tipole of I = 2000 (2-rrGp(z)a 2 c- 2 Wi( x)\/ Pss(k, x)) [b l ue do ttedl 
and the II integrand for the orig inal lffirata fc Seliakl l|2004l l LA 
model (Cip(z)arii (x)\/ Pss(k, x)) [ re d dotted] and the LA model 
given in this work (Cip(z = 0)n;(x)^/ Pgg n (k, x)) [red solid], both 
shown for I = 2000. We present integrands without [upper panel] 
and with the appropriate matter power spectrum term [lower 
panel]. Each angular power spectrum is plotted for the first and 
tenth tomographic redshift bins of our fiducial survey, these ap- 
pear to the left and right hand of each plot respectively. Here 
the term P(k) corresponds to either the linear or non-linear mat- 
ter power spectrum depending on the term (G or 7), and on the 
implementation used (HS04NL or This Work). 



creases with redshift in the HS04NL prescription. As IAs are 
normalised against a Ci value measured at low redshift this 
means that the new II term is comparable to the HS04NL 
term at low redshift (if it were possible to show the II terms 
at 2 = on the top panel of Fig. [4] they would be identical) 
but increasingly diverge a higher redshifts. This is why the 
II kernel for the 10th (high-z) redshift bin is reduced more 
as we move from HS04NL to the new model than the equiv- 
alent lines for the 1st (low-z) bin. As the shear-shear term 
is unaffected by our treatement of IAs this means that IAs 
make a less significant contribution to the total shear signal 
at higher redshift than in the HS04NL approach. 

The impact of moving from the common NLA approach 
to our motivated treatment of nonlinear clustering and IAs 
is to weaken the contribution of IAs relative to cosmic shear 
because we discard all nonlinear clustering power for the 
I A terms. This means we lose all power due to nonlinear 
clustering in the II term and half the nonlinear power from 
the GI term where the I contribution is linear but the G 
contribution remains dependent on the fully nonlinear mat- 
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Figure 5. Projected angular power spectra which contain IA contributions. Computed for our fiducial survey and displayed for a variety 
of tomographic redshift bin combinations as in Fig[T] Power spectra are computed using the HS04NL approach [narrow lines] and the 
new approach described in this work [thick lines] . The GG (solid black narrow) power spectrum is shown for comparison. 



ter power spectrum. The impact is felt at all redshifts but 
the effect is strongest at low-z where nonlinear clustering is 
strongest. Fig. [5] shows the projected angular power spectra 
for those components which are affected by the change from 
the HS04NL model (as used in Ijoachimi fc Bridle! j2009h ) 
to the latest implementation described in this work. The 
format of the plot is the same as Fig. [T] lines correspond- 



ing to the latest LA implementation are thick versions of 
the same line-style as their HS04NL equivalents. The shear- 
shear (GG) angular power spectra are, of course, unaffected 
by the change, they are shown in the upper triangle for com- 
parison. 
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